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UNSTABLE  WAVES  OF  JET  FLOWS  WITH 
DENSITY  INHOMOGENEITY 


i 


INTRODUCTION 

The  problem  to  be  considered  concerns  linear  unstable  waves  of  axisymmetric  jet 
flows  in  the  presence  of  density  inhomogeneities.  Such  flow  phenomena  occur  when  a  jet 
is  discharged  into  a  stratified  medium,  e.g.,  pollutants  and  industrial  waste  discharged 
into  the  environment,  cooling  water  discharged  from  power  plants  into  rivers  and  lakes, 
and  flow  patterns  generated  by  vehicles  moving  in  the  ocean.  Investigating  the  instability 
characteristics  of  such  jet  flows  is  necessary  to  fully  understand  the  overall  behavior  of  the 
flow  patterns  and  their  corresponding  effects  on  the  environment.  Changing  the  instability 
characteristics  of  the  flow  can  mean  controlling  the  flow  patterns,  for  example,  produced  by 
vehicles  moving  in  the  ocean.  As  in  the  case  of  two-dimensional  shear  flows,  axisymmetric 
jet  flows  may  possess  amplified  waves  due  to  the  Kelvin- Helmholtz  mechanism  except  that 
the  formulation  for  jet  flows  is  complicated  by  the  absence  of  the  Squire  transformation 
and  the  consideration  of  cylindrical  geometrv.  *  p£o  t-j/ 

In  a  well-known  paper  concerning  unbounded  parallel  flows  of  the  jet-wake  type. 
Batchelor  &  Gill  (1962)  presented  a  general  mathematical  analysis  of  the  stability  charac¬ 
teristics  of  axisymmetric  flows  of  homogeneous  fluids.  A  necessary  condition  for  instability 
and  a  semi-circle  theorem  for  possible  unstable  waves  were  derived.  By  discussing  the 
governing  equation  near  the  jet  axis,  they  were  able  to  reach  some  general  conclusions 
on  the  characteristic  difference  between  the  axisymmetric  and  nonaxisymmetric  modes. 
Meanwhile  Reynolds  (1962)  conducted  an  experimental  study  and  observed  axisymmetric 
condensations  and  puffs  in  a  water  tank.  The  axisymmetric  mode  of  a  jet  column  was  also 
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examined  by  Gill  (1962)  with  a  slightly  viscous  jet  flow  model  and  later  by  Mattingly  & 
Chang  (1974)  with  a  coordinated  theoretical  and  experimental  investigation.  The  latter 
found  that  the  dominant  disturbance  in  the  jet  was  an  axisymmetric  one.  It  is  expected 
that  initial  disturbances  will  follow  the  direction  of  mean  flows  and  the  axisymmetric  waves 
will  in  general  first  amplify  for  the  onset  of  instability. 

In  this  theoretical  study,  the  effect  of  radial  density  variations  on  the  stability  of 
incompressible  axisymmetric  jet  flows  will  be  investigated.  Both  the  necessary  condition 
for  instability  and  the  semi-circle  bound  on  amplified  waves  will  be  obtained  through  a 
unified  integral  representation.  The  results  obtained  here  show  that  the  necessary  condition 
for  the  existence  of  amplified  waves  depends  not  only  on  the  velocity  profile,  but  also  on 
the  effect  of  density  inhomogeneities.  Both  positive  and  negative  density  gradients  have 
stabilizing  effects.  The  semi-circle  theorem  for  unstable  waves  can  be  extended  to  flows 
with  radius-dependent  density.  For  top-hat  profiles  of  jet  or  wake  type,  the  semi-circle 
bound  is  found  to  be  the  best  possible. 


GOVERNING  EQUATIONS  AND  INSTABILITY  CHARACTERISTICS 


The  jet  flow  to  be  considered  consists  of  a  radius-dependent  axial  velocity  W(r)  in  a 
fluid  with  density  p(r),  and  is  confined  within  the  annular  region  between  two  concentric 
cylinders  located  at  r  =  Rx  and  r  =  R2.  Neglecting  the  gravitational  effects,  one  finds, 
within  the  framework  of  the  normal  mode  analysis,  that  sinusoidal  waves  are  governed  by 
the  two  ordinary  differential  equations 

2 

(k2  +  ~ )p  =  ipk[DW  -  (W  -  c)(Du  +  -)],  (1) 

r*  r 

Dp  =  —  ipk(W  —  c)u  (2) 

with  D  =  d/dr. 

Here  u  and  p  are  the  perturbations  in  radial  velocity  and  pressure  respectively.  To 
consider  only  temporal  instabilities,  the  azimuthal  wave  number  m  is  an  integer,  the  axial 
wave  number  k  is  real,  while  the  phase  velocity  c  =  cr  -I-  ic ,  is  in  general  complex.  Since 
the  equations  are  invariant  under  complex  conjugation,  a  non-zero  c,  implies  instability. 

Combining  equations  (1)  and  (2)  and  eliminating  the  variable  p.  one  finds  a  single 
differential  equation  governing  instability  as  follows: 

(W  -  c)^D[pq{Du  +  ^)j  -  puj  -  D.(pqDlV)u  =  0  to) 

where  D*  =  D  -  1/r  and  q  =  r2/(m2  4-  k2r2).  Equation  (3)  and  boundary  conditions 
u(Ri)  =  u(R2)  ~  0  represent  an  eigenvalue  problem  for  the  inviscid  unstable  waves  n 
axisymmetric  jet  flows  in  a  stratified  medium.  Reminiscent  of  the  Rayleigh  equation 
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encountered  in  the  stability  of  two-dimensional  shear  flow,  equation  (3)  possesses  a  first- 
order  singularity  for  neutral  disturbances  and  analytical  solutions  in  terms  of  well-known 
functions  are  in  general  not  possible  except  for  some  special  flow  profiles.  We  will  emphasize 
the  effect  of  density  inhomogeneities  and  seek  solutions  to  those  flow  profiles  which  admit 
exact  solutions  in  terms  of  modified  Bessel  Functions.  To  do  this,  one  must  avoid  the  first 
order  singularity  in  the  equation  as  c,  vanishes. 

Before  obtaining  exact  solutions  to  equation  (3)  we  will  briefly  re-examine  the  inviscid 
instability  characteristics  of  the  flow  via  an  integral  representation.  The  analysis  in  the 
present  case  is  particularly  simple  and  the  principal  steps  of  deriving  these  characteristics 
are  reproduced  here  for  the  sake  of  completeness.  To  perform  the  analysis,  we  look  for 
unstable  waves  with  c<  >  0  and  make  the  transformation  u  =  (W  —  c)”t/>.  Here  n  is  an 
arbitrary  real  constant.  Multiplying  the  resultant  equation  by  the  complex  conjugate  r<j> , 
integrating  it  over  the  flow  domain,  and  applying  the  boundary  conditions  that  <t>  vanishes 
at  the  solid  boundaries  of  the  inner  and  outer  cylinders,  one  obtains 

j*'{(W  -  c)2np[q\D*1>\2  +  M2]  -  (W  -  c)2n-l[(n  - \)D(pQ)}r\rP\2 

-{W  -  c)2"“2[n(n  -  l)p9(DW)2]M2  jrdr  =  0  (4) 

where  Q  «  rZ?W/(m2  +  k2r2)  and  D*  =  D  +  1/r.  Letting  n  —  0  one  finds,  from  the 
imaginary  part  of  the  resultant  integral 

In/d2 

«/„,  <5> 


For  instability,  the  expression  D(pQ)  must  vanish  at  least  at  one  interior  point  within  the 
flow  domain,  i.e., 


dq.dp 
Q  p 


=  o 


(6) 


The  first  term  of  the  above  equation  is  reminiscent  of  the  inflexion  point  theory  encoun¬ 
tered  in  plane  parallel  flows,  and  has  been  discussed  in  detail  in  the  study  of  homogeneous 
jet  flows  (Batchelor  &  Gill  1962).  The  second  term  depends  solely  on  density  inhomo¬ 
geneities  and  will  have  significant  contributions  when  the  condition  for  instability  of  the 
otherwise  homogeneous  flow  field  is  satisfied.  It  is  shown  from  equation  (6)  and  later  via 
the  exact  solution  for  a  special  flow  profile,  that  the  necessary  condition  for  amplified  waves 
established  for  homogeneous  fluids  is  not  sufficient  if  density  inhomogeneities  are  present. 
Cylindrical  vortex  sheets  can  be  stable  if  large  density  gradients  are  present  in  the  flow. 


The  refinement  of  condition  (6)  can  be  obtained  by  equation  (5)  and  the  real  parf  of 
equation  (4)  for  n  =  0.  One  can  easily  show  a  strong  necessary  condition  for  instability  is 


D{ 


prDW 
m 2  +  k2r2 


}{W  -  W,)  <  0 


>7) 


where  W,  is  the  mean  velocity  evaluated  at  a  certain  point  inside  the  flow  field  where 
condition  (6)  is  satisfied. 

The  semi-circle  bound  can  easily  be  derived  by  combining  the  real  and  imaginary  parts 
of  integral  (4)  for  n  =  1.  Imposing  the  upper  and  lower  limits  a  and  6  for  the  mean  axial 
velocity,  one  obtains 


/•** 

0>  /  (W  -  a)(W  -  b)p[q\Dm*\2  +  \4>\2}rdr 
JRi 

=  {[cr-  +  6)J2  +  c2-  [i(a  -  6))2}  p[q |DV|2  +  \rj>\ 2}rdr  (8a) 

[cr-^(a  +  6)]2+c?<[|(a-6)]2  (86) 

Thus  the  semi-circle  theorem  is  valid  in  the  presence  of  density  inhomogeneities,  saying 
that  the  complex  wave  velocity  must  lie  inside  a  semi-circle  with  the  diameter  equal  to 
the  range  of  the  velocity.  The  subsequent  section  will  show  that  this  semi-circle  bound 
provides  an  exact  solution  for  some  particular  flow  profiles. 


TWO  EXACT  SOLUTIONS  AND  THEIR  SEMI-CIRCLE  BOUNDS 


To  examine  the  validity  of  the  earlier  obtained  instability  characteristics  and  to  gain 
some  information  on  instability  growth  rates  as  a  function  of  velocity,  density  and  wave 
numbers,  one  must  seek  explicit  solutions  to  equation  (3)  for  some  spec: til  flow  configu¬ 
rations.  We  will  consider  those  profiles  for  which  solutions  in  terms  of  modified  Bessel 
functions  are  possible.  One  of  such  profiles  to  investigate  is 


w 

=  W,, 

for 

R\  <r  <  R 

w 

£ 

ii 

for 

R<r  <  R? 

where  W\,  Wi,pi,p2,(Ti  and  cti  sue  constants.  This  flow  represents  a  top-hat  jet  core  with 
radius-dependent  density  surrounded  by  a  fluid  column  with  different  density.  Because  the 
discontinuity  in  the  velocity  at  the  interface  r  =  R  satisfies  instability  condition  (6)  for 
homogeneous  fluids,  this  velocity  profile  will  generally  be  unstable  for  all  modes  in  uniform 
media  via  the  Kelvin- Helmholtz  mechanism.  However,  it  will  be  demonstrated  later  that 
these  unstable  waves  can  be  stabilized  if  large  density  gradients  are  considered. 

The  perturbations  in  radial  velocity  and  pressure  for  the  flow  profile  in  equation  (O' 

are  given  by 


u,  =  r 
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I„t{kr) 
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(10) 


Pj  =  -ipjR-°>  k(W}  -  c)r'i{AjIVj(kr)  +  B}Ku>{kr)} 


with  uj  =  y/m*  +  (<Tj/ 2)2  for  j  =  1,2.  Here  J„,(*r)  and  Ku>(kr)  are  the  modified  Bessel 
functions  of  the  first  and  second  kind  of  order  v},  and  a  prime  denotes  the  total  derivative 
with  respect  to  the  quantity  shown. 

By  integrating  equations  (1)  and  (2)  across  the  interface  at  r  =  R ,  both  the  kinematical 
and  dynamical  interfacial  conditions  are  obtained  as  follows: 


>=  0 


where 


fy-c 

<  p  >=  0 

<  4>  >=  —  (H_o). 
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Making  use  of  the  boundary  conditions  at  r  =  R\ ,  R2  and  applying  the  matching  conditions 
(11)  yield  the  secular  relation  for  instability 

(Wi  -  c)[IVi(k)  -  HiK^jK)) 


—a 


(W2  -c )[#„(«)  -  H2IVi(K)\ 


r(72  ,  *-&[,,(*),  r,  tr  ,  K'W'C 


-  =  0 


(12) 


-}K„7{k)  -  4- 


Here 
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|f|/2  (  ) 


h2  = 


°  A„1(/c2) _ 


rCT2  .  K2Iu,(K2)lr  ,  , 

hr  +  Tfr^1''2' 
1  u,(*2) 


where  a  =  p2/pi,  k  =  fcH,  =  fcHj,  and  «2  =  fcH2.  The  expression  Hx  vanishes  as  Rx  —  0 
and  so  does  H2  as  R2  —  oo  to  represent  solutions  for  unbounded  flows.  Solving  equation 
(12)  for  the  wave  velocity  yields 


c  = 


Ex  Wx  -  aE2  W2  £  ( Wj  -  \V2 )( aE\  Ei )  * 
E\  —  afi2 


Here 


Ei  = 


[IM-HiK^k)] 
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It  can  be  shown  that  E\  >  0  and  E?  <  0  for  all  the  values  of  «i  <  x  <  *c2-  Nonzero  c, 
is  therefore  expected  for  all  inodes  except  for  very  small  or  large  density  ratios.  A  liquid 
jet  impinging  into  the  atmosphere  is  of  course  stable  for  all  modes  even  in  the  presence 
of  the  strong  Kelvin-Helmholtz  effects  arising  from  the  cylindrical  vortex  sheets.  This 
stabilization  phenomenon  stems  from  the  fact  that  instability  condition  (6)  is  violated  by 
large  density  gradients  even  though  DQ  vanishes  inside  the  flow  field.  One  other  interesting 
feature  we  like  to  point  out  is  that  all  amplified  waves  marked  by  equation  (13),  in  spite  of 
the  density  variations,  wave  numbers  and  locations  of  the  solid  boundaries,  lie  exactly  on 
the  boundary  described  by  equation  (8b)  with  the  stable  solutions  corresponding  to  a  =  0 
and  a  =  oo  at  both  ends  of  the  semi-circle.  This  exactness  may  be  explained  by  examining 
integral  (8a)  for  the  flow  under  consideration.  The  velocity  profile  for  the  two-region  flow 
can  be  written  as 


W  =  Wi  +  (W2  -  Wi  )H(r  -  R)  for  Ri  <  r  <  R2  (14) 

where  H(r— R)  is  the  Heaviside  function.  The  equality  of  equation  (8b)  is  obtained  through 
the  substitution  of  equation  (14)  for  W  in  integral  (8a).  Note  that  the  instability  growth 
rate  for  fixed  k  first  increases  with  increasing  a,  reaches  a  maximum  at  an  intermediate 
value  and  then  decreases  with  still  larger  a,  with  the  maximum  growth  rate  (|W2  —  Wi|/2) 
governed  by 
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(15) 


The  special  case  of  01  =  02  =  0  in  equation  (15)  for  unbounded  flows  is  plotted  in  figure 
(1)  which  shows  the  characteristic  difference  between  the  axisymmetric  mode  and  the 
asymmetric  ones. 


The  second  example  whose  instability  characteristics  we  wish  to  study  for  axisymmet¬ 
ric  perturbations  is  a  three-region  flow  with  the  profile  given  as: 
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P-  P3 

R2  <  r  <  co 

Here  pi,p2,ps  and  a  are  arbitrary  constants.  This  profile  can  be  used  to  model  a  jet 
exhausting  into  an  environment  of  different  density  and  velocity  with  the  middle  region 
representing  the  shear  layer  to  allow  transitions  of  the  physical  quantities  of  the  jet  to 
its  surroundings.  One  can  examine  such  flow  instabilities  by  adjusting  the  parameters 
Wj,Rj,pj  and  <r.  The  perturbation  radial  velocities  and  pressures  for  the  flow  in  the  inner 
and  outer  regions  can  be  reduced  from  equation  (10)  by  setting  m  =  0  and  =  0.  while 
the  solutions  in  the  middle  region  tire  given  by 

U2  =  4/M(fcr)  4-  BK„{kr)} 


P2  =  $  )r2‘  W(*>->  +  BK.itr)] 

-l$^Ww'+W^w’-c 

(IT) 

where  p  =  1  —  <r/2.  Applying  the  matching  conditions  (11)  at  both  interfaces  r  =  /?,  and 
r  as  R2  gives  the  secular  relation 
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Rearranging  determinant  (18)  by  using  the  identity 

WKll(2)-i*(2)Ki(z)  =  K 

one  finds  the  solution  for  the  wave  velocity 

_Wl  +  W2  ii(W2-W1)(F1+Gi6*»)Ill(Kl)K»(K2)-(Gi  +  F26^)K^Kl)Itl(K2) 
r  2  1  -6**  F1G2Im(k1)Kii(k2)-G1F2Kii(k1)Iii(k2) 


=  fx(W2-W,) _ n/A _ 

1  —  62*  F\G2Iu(k\  )Kii(k2)  —  GiF2Klt(Ki)Ili(K2) 

A  =  462"  -  ( i— +  -±-  -  V m  CW*!  )K»(k2  ) 


(19a) 

(196) 


1 

1 

1 

s2* 

2p 

+  c2 

"K 

1 

1 

i 

62tt 

2  H 

f2 

~  g7 

“ ( — 2^ —  +  ^  -  ^)Gi^(*i)4(K2)J  (19c) 

Thus  amplified  waves  exist  whenever  A  >  0.  equations(19)  for  the  particular  case  W2  = 
0,  Pi  =  fh,  —  P3  and  <7  =  0  reduce  to  the  result  investigated  earlier  by  Michalke  &:  Schade 
(1963)  in  their  study  of  shear  flow  instability  in  uniform  fluids  except  that  0*  should  be 
dropped  from  their  equations  (91)  and  (94).  It  should  be  pointed  out  that  the  complex 
wave  speed  in  their  limiting  case  as  Ri  approaches  R2  appears  to  be  incorrect.  The  solution 
for  unbounded  cylindrical  vortex  sheets  should  have  been  recovered  if  the  proper  limiting 
process  had  been  taken.  The  flow  in  that  limiting  case  was  certainly  unstable  for  all  axial 
wave  numbers  with  the  real  wave  velocity  approaching  the  velocity  of  the  center  of  the  jet 
as  k  — ►  0.  The  limiting  process  as  jRi  — *  R2  reduces  equations(19)  to 


W!+W2  (W2  -  W2)  [/„(*)#!(*)  -  aAV*)/,^)] 

2  2  [I0(k)Ki(k)  +  oiKo(k)Ii(k)} 


(20  a) 


c,  =  ±(W,  -  W2) 


\J  qIq(k)Kq(k)/i(k)Ki(k) 
I0{k)Ki{k)  +  aA'0(K)/i(/c) 


with  a  representing  the  density  ratio  between  the  outer  and  inner  regions,  equations  (20) 
can  of  course  be  reduced  from  equation  (13)  for  the  special  case  of  o\  =  <r2  =  0  and  m  =  0 
for  unbounded  axisymmetric  waves.  It  in  fact  can  be  shown  that  the  general  expression 
(19a)  for  cr  possesses  a  limiting  value  W2  for  long  axial  wave  lengths.  This  acts  as  a 
support  to  Batchelor  &  Gill’s  finding  that  long  waves  with  axisymmetry  travel  with  the 
speed  of  the  center  of  the  jet.  The  corresponding  instability  discriminant  reduced  from 
(19c)  takes  the  form 

A  = - - — r[oip(l  -  62(t)  +  (1  -  62,t)2}K2logK2  >  0  (21) 

Ot\Ct2H* 
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Thus  long  waves  are  always  unstable  except  for  k  =  0.  For  sufficiently  large  wave  numbers, 
the  asymptotic  form  of  the  modified  Bessel  functions  reduces  equation  (19c)  to 


A  =  -[gijk±Mgipfk£^  <  0  (22) 

Therefore  amplified  waves  cannot  exist  in  their  short  wavelength  range  except  for  the 
limiting  case  f?i  =  R2  in  which  the  maximum  growth  rates  for  uniform  fluids  occurs  in 
the  limit  k  — ►  00.  These  instability  behaviors  axe  plotted  out  in  Figs.  (2)  and  (3)  for  the 
special  case  of  aj  =  a?  =  <7  =  1  and  W2  =  0,  with  «  =  R2  —  R\  representing  the  shear  layer 
thickness.  Both  the  velocity  and  density  in  this  case  vary  linearly  across  the  shear  layer. 
Figure  (2)  shows  the  instability  growth  rates  in  comparison  with  the  eigenvalue  bound  in 
their  complex  velocity  domain.  The  amplified  rate  for  fixed  shear  layer  thickness  e  first 
increases  with  increasing  k,  reaches  a  maximum  value  at  an  intermediate  value  of  wave 
number,  and  decreases  to  zero  at  the  relatively  large  value  of  k  given  by  equation  (19c) 
when  the  discriminant  vanishes.  These  instability  growth  rates  approach  the  semi-circle 
bound  as  the  shear  layer  thickness  diminishes  and  the  maximum  growth  rate  \W2  —  W\\/2 
for  cylindrical  vortex  sheets  is  reached  when  e  =  0.  Note  that  the  complex  wave  velocity 
expands  to  the  left  portion  inside  the  semi-cir  ;lt  as  the  density  ratio  between  the  jet 
core  and  the  surroundings  decreases.  For  sufficiently  small  or  large  density  ratios,  the 
flow  is  stable  against  all  disturbances.  The  neutral  stability  boundary  and  the  maximum 
growth  rates  are  shown  in  Figure  (3)  as  a  function  of  the  wave  number  and  the  shear  layer 
thickness.  The  figure  suggests  that,  for  non- zero  shear  layer  thickness,  disturbances  for 
which  e  >  2/k  are  stable,  and  the  maximum  growth  rate  occurs  when  ke  is  of  the  order  of 
one-tenth.  It  is  also  shown  through  this  particular  flow  configuration  that  amplified  waves 
with  axisymmetry  may  exist  for  some  slowly  varying  profiles  in  the  long  axial  wavelength 
range. 


CONCLUDING  REMARKS 


The  instability  characteristics  of  axisymmetric  jet  flows  in  inhomogeneous  fluids  were 
discussed.  The  exact  solutions  to  the  governing  stability  equation  confirm  these  char¬ 
acteristics;  in  particular,  the  semi-circle  theorem  provides  the  best  possible  bound  on  all 
amplified  waves  of  the  top-hat  type  velocity  profiles.  The  discharge  of  a  jet  into  a  stratified 
medium  may  produce  some  organized  flow  patterns  which  do  not  exist  in  a  homogeneous 
environment.  An  understanding  of  the  instability  characteristics  is  of  practical  interest, 
because  changing  those  characteristics  of  the  flow  can  mean  controlling  the  flow  patterns, 
for  example,  produced  by  vehicles  moving  in  the  ocean. 
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Figure  1.  Curves  of  Maximum  Growth  Rates  <\IV2  -  IF,  ,72)  for 


'7\ 
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elocity 


Figure  3.  Variations  of  the  Stability  Boundary  and  the  Maximum  Growth  Rates  for  the  Special 
Case  oi  =  a2  =  <7  =  1  and  Wi  =  0. 
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